EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2012-361 
2012/12/27 



CMS-SUS-12-017 



Search for new physics in events with same-sign dileptons 
and b jets in pp collisions at y/s = 8 TeV 



The CMS Collaboration 



Abstract 



A search for new physics is performed using events with isolated same-sign leptons 
and at least two bottom-quark jets in the final state. Results are based on a sample of 
proton-proton collisions collected at a center-of-mass energy of 8 TeV with the CMS 
detector and corresponding to an integrated luminosity of 10.5 fb _1 . No excess above 
the standard model background is observed. Upper limits are set on the number of 
events from non-standard-model sources and are used to constrain a number of new 
physics models. Information on acceptance and efficiencies is also provided so that 
the results can be used to confront an even broader class of new physics models. 
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1 Introduction 



Events with same-sign, isolated leptons (e or y) from standard model (SM) processes in proton- 
proton collisions are extremely rare. As a result, searches for anomalous production of same- 
sign dileptons can be very sensitive to new physics contributions. These include supersymme- 
try (SUSY) HH3, universal extra dimensions [4j, pair production of T5/3 particles (fermionic 
partners of the top quark) [5], heavy Majorana neutrinos HI, and same-sign top-pair produc- 
tion HI. More specifically, the signature of same-sign dileptons accompanied by bottom 
quarks can arise from SUSY processes where third-generation quark superpartners are lighter 
than other squarks |(9 HT3|I , resulting in an abundance of top and bottom quarks produced in cas- 
cade decays. These scenarios are particularly interesting for two reasons. First, because they 
have not been excluded by early SUSY searches at the Large Hadron Collider (LHC); second, 
because relatively light bottom and top squarks are necessary if SUSY is to be the "natural", 
(i.e., not fine-tuned), solution to the gauge hierarchy problem (see for example the discussion 



In this Paper we present the result of a search for new physics in events with same-sign iso- 
lated dileptons and two or more bottom-quark jets (b jets). The search was performed on a 
dataset corresponding to an integrated luminosity of 10.5 fb _1 of proton-proton (pp) collisions 
collected at a center-of-mass energy of 8 TeV with the Compact Muon Solenoid (CMS) [14] de- 
tector at the LHC. This search is an extension of the search at 7 TeV described in Ref. [15|. 

In the analysis, we select events with two isolated, high transverse momentum (pj), same- 
sign leptons (e or y), and at least two jets identified as b jets (b-tagged jets). The numbers 
of observed events in signal regions defined by additional requirements on missing transverse 
energy (E™ 1SS ), the number of b-tagged jets, and the scalar sum of jet pj (Hj), are compared with 
SM expectations. Having found no evidence for a beyond-the-SM contribution to the event 
counts, we set limits on the number of events from new physics in the event sample. These 
limits are then translated into a limit on the same-sign top-pair production cross section and 
are used to bound the parameters of other new physics models. Finally, we include additional 
information on the event selection efficiencies so that our results can be used to test models of 
new physics not explicitly considered here. 



The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. The CMS experiment uses a right-handed coordinate 
system, with the origin defined to be the nominal interaction point, the X axis pointing to the 
center of the LHC ring, the y axis pointing up (perpendicular to the LHC plane), and the z 
axis pointing in the counterclockwise beam direction. The polar angle 9 is measured from the 
positive z axis and the azimuthal angle (p is measured in the x-y (transverse) plane. The pseu- 
dorapidity t] is defined as tj = — In [tan (8/2)}. Within the field volume are a silicon pixel and 
strip tracker, a crystal electromagnetic calorimeter, and a brass-scintillator hadron calorime- 
ter. Muons are measured in gas-ionization detectors embedded in the steel return yoke. Full 
coverage is provided by the tracker, calorimeters, and the muon detectors within \rj\ < 2.4. 
In addition to the barrel and endcap calorimeters up to | rj | =3, CMS has extensive forward 
calorimetry reaching \tj\ < 5. A more detailed description of the CMS apparatus can be found 
in Ref. El. 



in Ref. (111). 



2 The CMS detector 
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4 Backgrounds 



3 Event selection and Monte Carlo simulation 

The event selection is essentially the same as for the 7TeV analysis |15[. We require two isolated 
same-sign leptons (e or y) with pair invariant mass above 8 GeV, and at least two b-tagged 
jets. Events are collected with dilepton triggers. The kinematic requirements on leptons and 
jets are summarized in Table [Tj Events with a third lepton are rejected if the lepton forms an 
opposite-sign same-flavor pair with one of the first two leptons, and if the invariant mass trigg 
of the pair satisfies w% < 12GeV (pj > 5GeV) or 76 < < 106 GeV (pj > 10 GeV), where 
the conditions in parentheses show the respective pi requirement on the third lepton. These 
requirements are designed to minimize backgrounds from events with 7* — > £ + £~ , low-mass 
bound-state, and multi-boson (WZ and ZZ, as well as tribosons) production. 



Table 1: Kinematic requirements for leptons and jets in this analysis. 





Pt 


n 


electrons 

muons 

jets 

b-tagged jets 


p T > 20 GeV 
p T > 20 GeV 
p T > 40 GeV 
p T > 40 GeV 


\r]\ < 1.442 or 1.566 <\rj\< 2.4 
\t}\ < 2.4 
\t]\ < 2.4 

\v\ < 2A 



With respect to Ref. [15], there have been a few modifications in the lepton identification and 
isolation requirements; the main difference is that lepton isolation is now computed with 
particle-flow information [16] rather than with charged tracks and calorimeter deposits, and 
an event-by-event correction is made to account for the effect of multiple pp interactions in 
the same bunch crossing. This correction is based on a measurement of the hadronic event 
activity not associated with the pp interaction that produced the leptons (the mean number of 
interactions in the same crossing is 15). 

The jet reconstruction is the same as in Ref. |[T5l , but the b-tagging algorithm has been changed 
slightly. We now use the CSV method of Ref. [17], which is based on the combination of 
secondary-vertex reconstruction and track-based lifetime information. Efficiencies for lepton 
identification, trigger efficiency, and b-tagging are summarized in Sections [6] and [7[ 

Monte Carlo (MC) simulations are used to estimate some of the SM backgrounds (see Sec- 
tion El, as well as to calculate the efficiency for various new physics scenarios. In this analysis, 
SM background samples and new physics signal samples are generated with the MadGraph 
V5 fl8j and PYTHIA V6.4 [19] programs, respectively. Standard model processes are simulated 
using a GEANT4-based model [20] of the CMS detector; the simulation of new physics signals 
is performed using the CMS fast simulation package (2l|. All simulated events are processed 
with the same chain of reconstruction programs used for collision data. 

4 Backgrounds 

There are three sources of SM background to potential new physics signals: 

• "fake leptons", i.e., leptons from heavy-flavor decay, misidentified hadrons, muons 
from light-meson decay in flight, or electrons from unidentified photon conversions. 
These are estimated from a sample of events with at least one lepton that passes a 
loose selection but fails the full set of tight identification and isolation requirements. 
The background estimation uses the "tight-to-loose" ratio, i.e., the probability for a 
loosely identified fake lepton to also pass the full set of requirements, as determined 
from studies of fake leptons in multijet events. 
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• "charge flips", i.e., events with opposite-sign isolated leptons where one of the lep- 
tons is an electron and its charge is misreconstructed due to severe bremsstrahlung 
in the tracker material (this effect is negligible for muons). Charge flips are esti- 
mated by selecting opposite-sign ee or eji events passing the full kinematic selection, 
weighted by the pj- and //-dependent probability of electron charge misassignment. 
This probability is obtained from MC simulation and is validated on a sample of 
Z — > ee events. 

• rare SM processes that yield same-sign high-px leptons and b jets, mostly from tfW 
and tfZ production. This background contribution is obtained from Monte Carlo 
simulation. 

More details on the background estimation can be found in Refs. fl5H22fl . The systematic un- 
certainties are 20% for charge flips, and 50% for both fake and rare SM backgrounds. The next- 
to-leading order ttW and ttZ production cross sections used to normalize the MC predictions 
are 232 fb ESI and 208 fb HM&J, respectively. 

5 Event yields 

The search is based on comparing the event yields and background predictions in nine signal 
regions (SRs) with different requirements on E™ ISS , Hj, and the number of b-tagged jets. In 
order to improve the sensitivity to pp — > tt production, we also define a SR with positive 
leptons only. The requirements for the SRs are loosely motivated by possible new physics 
scenarios. For example, signal regions with high EJp lss are most sensitive to R-parity-conserving 
SUSY models with relatively light neutralinos; requiring larger Hj reduces the SM background 
while generally maintaining good efficiency to detect new heavy particles that are produced 
via the strong interaction; the signal region with three or more b-tagged jets is sensitive to final 
states that include several b quarks. 

Table 2: A summary of the results of this search. For each signal region (SR), we show its 
distinguishing kinematic requirements, the prediction for the three background (BG) compo- 
nents, as well as their total, and the observed number of events. We also show the 95% CL 
upper limit on the number of new physics events (Nm) under different assumptions on the 
signal efficiency uncertainty. Note that the number of jets on the first line of the table includes 
both b-tagged and non b-tagged jets. 





SRO 


SRI 


SR2 


SR3 


SR4 


SR5 


SR6 


SR7 


SR8 


No. of jets 


> 2 


> 2 


> 2 


> 4 


> 4 


> 4 


> 4 


> 3 


> 4 


No. of btags 


> 2 


> 2 


> 2 


> 2 


> 2 


> 2 


> 2 


> 3 


> 2 


Lepton charges 


+ + /-- 


+ + /-- 


±± 


+ + /-- 


+ + /-- 


+ + /-- 


+ + /-- 


+ + /-- 


+ + /-- 


£miss 


>0GeV 


>30GeV 


>30GeV 


>120GeV 


>50GeV 


>50GeV 


>120GeV 


>50GeV 


>0GeV 


Ht 


>80GeV 


>80GeV 


>80GeV 


>200GeV 


>200GeV 


>320GeV 


>320GeV 


>200GeV 


>320GeV 


Fake BG 


25 ±13 


19 ±10 


9.6 ±5.0 


0.99 ± 0.69 


4.5 ± 2.9 


2.9 ±1.7 


0.7 ±0.5 


0.71 ± 0.47 


4.4 ± 2.6 


Charge-flip BG 


3.4 ±0.7 


2.7 ±0.5 


1.4 ±0.3 


0.04 ±0.01 


0.21 ± 0.05 


0.14 ± 0.03 


0.04 ± 0.01 


0.03 ± 0.01 


0.21 ± 0.05 


Rare SM BG 


11.8 ±5.9 


10.5 ±5.3 


6.7 ± 3.4 


1.2 ±0.7 


3.4 ± 1.8 


2.7 ±1.5 


1.0 ±0.6 


0.44 ± 0.39 


3.5 ±1.9 


Total BG 


40 ±14 


32 ±11 


17.7 ±6.1 


2.2 ±1.0 


8.1 ±3.4 


5.7 ±2.4 


1.7 ±0.7 


1.2 ±0.6 


8.1 ±3.3 


Event yield 


43 


38 


14 


1 


10 


7 


1 


1 


9 


Nm (13% unc.) 


27.2 


26.0 


9.9 


3.6 


10.8 


8.6 


3.6 


3.7 


9.6 


N UL (20% unc.) 


28.2 


27.2 


10.2 


3.6 


11.2 


8.9 


3.7 


3.8 


9.9 


Nm (30% unc.) 


30.4 


29.6 


10.7 


3.8 


12.0 


9.6 


3.9 


4.0 


10.5 



The definition of the signal regions and the observations of the search are summarized in Ta- 
ble [2] The basic kinematic properties of the selected events are shown in Fig. [Tj The results are 
consistent with the background predictions. 

Table [2] also shows the 95% confidence level (CL) observed upper limit (Nm) on the number 
of non-SM events, calculated using the CL S method |26l - l2"8| under three different assumptions 
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6 Efficiencies and associated uncertainties 
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Figure 1: Top plot: distribution of E™ 1SS versus Hj for the 43 events in the most inclusive signal 
region (SR0). Note that the > 2 jets requirement in SR0 implies Hj > 80 GeV. Bottom-left plot: 
projection of the top plot on the Hj axis. Bottom-right plot: projection of the top plot on the 
E™ 1SS axis. For the one-dimensional distributions, the number of events in each bin is scaled 
appropriately to reflect units of events per 10 GeV and is compared with the background pre- 
dictions, with their uncertainties. The binning of the one-dimensional distributions is chosen 
to match the E™ 1SS and Hj boundaries of the different signal regions. 



for the signal efficiency uncertainty. As can be seen from Table |2j the limits are only weakly 
dependent on this uncertainty. The uncertainty on the signal efficiency is discussed in Section|6] 

6 Efficiencies and associated uncertainties 

The dilepton triggers require one lepton to have pj > 17 GeV and the other pj > 8 GeV. The 
trigger efficiency, measured using events collected with hadronic triggers, is 96 ± 5% for ee, 
93 ± 5% for eji, and 88 ± 4% for jiji. All acceptances, which are calculated from MC simulation, 
include a correction factor for the trigger efficiency. 

The lepton selection efficiencies, determined using a sample of simulated events from a typical 
signal scenario involving multiple top quarks in the final states arising from gluino decay fl29l 
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CMS Simulation, is = 8 TeV CMS Simulation, vs = 8 TeV 




P T (GeV) p T (GeV) 

Figure 2: Lepton selection efficiency as a function of pi (left); b-tagging efficiency as a function 
of the b-quark pi (right). The superimposed curves represent the efficiency parametrizations 
discussed in Section [7] 

l30l , are displayed in Fig. [2] Studies of large samples of Z — > ee and Z jiji events indicate that 
the lepton selection efficiencies are reproduced to better than 2% by the simulation. However, 
leptons in SUSY events have significantly lower isolation efficiency due to the large amount of 
hadronic activity that is expected in these scenarios. As a result, we conservatively assign an 
uncertainty of 5% per lepton to the lepton selection efficiency. 

Figure |2] also shows the b-tagging efficiency, obtained from simulation, for b quarks with \f]\ < 
2.4. The uncertainty on this efficiency is a function of b-jet pi; it is 4.4% at pi = 40 GeV, 2.7% at 
pi = 100 GeV, and 10% at p T = 500 GeV B7I . 

The uncertainty on the energy scale of jets in this analysis is also a function of pi [31 J. It varies 
between 5% and 2% in the pi range 40-100 GeV for jets with \rj\ < 2.4. 

The importance of these effects depends on the signal region and the model of new physics. In 
general, models with high hadronic activity and large E™ ss are less affected by the uncertainty 
on the jet energy scale. 

The total uncertainty on the acceptance for the models described in Section [8] is in the 13-20% 
range. Finally, there is a 4.4% uncertainty on the yield of events from any new physics model 
due to the uncertainty in the luminosity normalization II32II . 



7 Information for model testing 

Our results can be used to confront models of new physics in an approximate way through 
generator-level studies that compare the expected numbers of events with the upper limits 
from Table |2j The prescription to be used is given in Ref. IITBl , Section 7. The E2p iss and Hi 
turn-on curves in this analysis are the same as those of Ref. [15]. However the lepton and b- 
tagging efficiencies are slightly different because of the changes in the underlying selections. 
Their current parametrizations are given below. 

The parametrization of the lepton efficiency, shown in Fig. [2] (left), is 



e = £00 erf 



p T -20GeV 



£20 



erf 



pi -20 GeV 



(1) 
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8 Limits on models of new physics 




Figure 3: Diagrams for the four SUSY models considered (Al, A2, Bl, and B2). 



with eoo = 0.65 (0.69), e 20 = 0.35 (0.48), and a = 42GeV (25GeV) for electrons (muons). 
The parametrization of the simulated b-tagging efficiency, shown in Fig. [2] (right), is e = 0.71 
for 90 < px < 170 GeV; at higher (lower) pj it decreases linearly with a slope of —0.0004 
(-0.0047) GeV" 1 . 

8 Limits on models of new physics 

The results of the search are used to constrain specific models of new physics. For each model 
considered, limits are derived from the signal region expected to give the most stringent limit 
on the cross section at a given point in the parameter space of the model. The event selection 
efficiency for a given model is obtained from MC simulation, and the limits are calculated 
including systematic uncertainties on lepton efficiency, trigger efficiency, luminosity, jet energy 
scale, and b-tagging efficiency. The latter two uncertainties are evaluated at each point in the 
parameter space. 

The results from SRI and SR2 are used to set limits on the cross section for same-sign top-quark 
pair production, 0"(pp — > tt + tt) from SRI, and c(pp — > tt) from SR2. Here c"(pp — > tt + tt) 
is shorthand for the sum c(pp — > tt) + 0"(pp — >■ tt). Note that in most new physics scenarios 
pp — > tt is suppressed with respect to pp — > tt because of the parton distribution functions of 
the proton. These limits are calculated using simulated pp — > tt events to model the acceptance. 
This acceptance, including branching fractions, is 0.29 ± 0.04%. We find upper limits c(pp — y 
tt + tt) < 0.87 pb and a(pp -> tt) < 0.30 pb at the 95% CL; the median expected limits are 0.72 
and 0.37 pb, respectively. 

Next, we present limits on the parameter spaces of four R-parity-conserving SUSY models 
with third-generation squarks. The decay chains under consideration are shown schematically 

in Fig- H 

Scenarios Al and A2 represent models of gluino pair production resulting in the tttt^j 1 ^ final 
state, where is the lightest neutralino [l2l|29H30H33H34l. m model Al, the gluino undergoes a 
three-body decay g — > tt^J mediated by an off-shell top squark. In model A2, the gluino decays 
to a top quark and and an anti-top squark, with the on-shell anti-squark further decaying into 
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an anti-top quark and a neutralino. 

Models Bl and B2 have bottom squarks decaying as bi — > tx,\ an d X\ ~ > W Model Bl 
is a model of bottom-squark pair production, followed by one of the most likely decay modes 
of the bottom squark; model B2 consists of gluino pair production followed by g — > bib. The 
gluino decay modes in models Al and A2 would be dominant if the top squark is the lightest 
supersymmetric quark. Conversely if the bottom squark is lightest, the decay mode in model 
B2 would be the most probable. 

Excluded regions in the parameter space of the four models are shown in Fig. |4| For the gluino- 
initiated models (Al, A2, and B2), we exclude gluinos with masses below around 1 TeV, with 
only a small dependence on the details of the models. This is because the limits are driven by 
the common gluino pair production cross section. In the case of bottom squark pair production 
(model Bl) we set a limit m(b t ) > 450 GeV at 95% CL. 

These models are also probed by other CMS new physics searches in different decay modes, 
although other searches have so far been interpreted in the context of model Al but not A2, 
Bl, or B2. For model Al the limits given here are complementary to the limits from the single- 
lepton and "razor" searches [36 J: less stringent at low m(x^) but more stringent at high m(J^{)- 
A similar conclusion applies to model A2, since the final state is the same. Comparable limits 
for model Al, as well as for similar models with top and bottom quarks from gluino decays, 
have been reported by the ATLAS collaboration H371439X In the case of bottom-squark pair 
production, the ATLAS collaboration has reported a limit m (bi ) > 390 GeV, but assuming the 
decay mode bi — > instead of the model Bl mode bi — > t%i K3- 

9 Summary 

In summary, we have presented results of a search for same-sign dileptons with bottom-quark 
jets using the CMS detector at the LHC, based on a 10.5 fb _1 data sample of pp collisions at ^/s 
= 8 TeV. No significant deviations from the standard model expectations are observed. 

The results are used to set upper limits on the same-sign top-pair production cross section 
c-(pp tt + tt) < 0.87 pb and cr(pp -> tt) < 0.30 pb at the 95% confidence level. 

We also exclude gluinos with masses up to approximately 1 TeV if they decay exclusively into 
top or bottom squarks, and we place a lower limit on the bottom-squark mass of 450 GeV. Our 
results assume that the top and bottom squarks decay as ti — > t^5 and bi — > tX\> respectively. 
In the latter case we have also assumed %7 ~~ ^ ^ %\- 
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Figure 4: Exclusion regions at 95% CL in the planes of m(^) vs. m(g) (model Al), m(x 1 ) 
vs. m(bi) (model Bl), m(ti) vs. m(g) (model A2), and m(bi) vs. m(g) (model B2). Models 
A2, Bl, and B2 have more than two mass parameters, and cannot be fully represented in a 
two dimensional plot. The assumed values of the additional mass parameters are indicated 
in the plots. The black lines represent the kinematic boundaries of the models. The excluded 
regions are those within the kinematic boundaries and to the left of the bands. The effects of 
the theoretical uncertainties on the next-to-leading-order plus next-to-leading-log calculations 
of the production cross sections f35l are indicated by the shaded bands; the expected limits and 
their ±1 standard-deviation variations are also shown. 
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